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ABSTRACT 

Quantum gravity effects are lilcely to play a crucial role in determining the outcome of 
gravitational collapse during its final stages. In this contribution we will outline a canonical 
quantization of the LeMaitre-Tolman-Bondi models, which describe the collapse of spher- 
ical, inhomogeneous, non-rotating dust. Although there are many models of gravitational 
collapse, this particular class of models stands out for its simplicity and the fact that both 
black holes and naked singularity end states may be realized on the classical level, depend- 
ing on the initial conditions. We will obtain the appropriate Wheeler-DeWitt equation and 
then solve it exactly, after regularization on a spatial lattice. The solutions describe Hawk- 
ing radiation and provide an elegant microcanonical description of black hole entropy, but 
they raise other questions, most importantly concerning the nature of gravity's fundamental 
degrees of freedom. 
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I. INTRODUCTION 



It is a great pleasure to have been invited to contribute to a festschrift in honor of Joshua 
N. Goldberg. His gentle manner, his fine intelligence and his modest nature are recognized 
by all. They conceal his ardor in defending the moral ideals he possesses and loves. 

So long as no generally agreed upon theory of quantum gravity exists, it is important 
to examine the quantization of particular models. This contribution addresses the quanti- 
zation of the LeMaitre-Tolman-Bondi (LTB) solutions, which describe the classical collapse 
of spherically symmetric, inhomogeneous, self-gravitating dust. The classical solutions were 
originally introduced by G. LeMaitre to study cosmology, where it has found interesting 
applications 0. Our principal interest here will be to develop the Hamiltonian formalism 
for both the classical and quantum LTB models, as such we will describe a generalization 
of work by Kuchaf P, who developed a midisuperspace quantization of the Schwarzschild 
black hole. While it would be preferable to take a fundamental field (e.g. a scalar field) 
for the matter part, this would make the formalism much less tractable Moreover, the 
relevant features of gravitational collapse already exhibit themselves for the dust model in 
the sense that the dust collapse may result classically in the formation of a black hole or of 
a naked singularity, depending on the initial conditions. 

In section II we review classical LTB collapse and present the canonical formalism for these 
models. The hypersurface action yields two constraints, viz., the Hamiltonian constraint 
and the momentum constraint. We reconstruct the mass and time from the canonical data 
and this leads naturally to new variables viz., the mass function, the dust proper time, the 
physical radius and their conjugate momenta, which are introduced along with the generator 
of the canonical transformation from the old to the new variables. The momentum conjugate 
to the mass function may be eliminated in the Hamiltonian constraint using the momentum 
constraint. This leads to a new and simpler constraint that is able to take the place of 
the original Hamiltonian constraint. We apply Dirac's quantization program to the new 
constraints in section HI and obtain the Wheeler-DeWitt equation. We then introduce a 
lattice regularization of the functional equations and find exact solutions, which we use in 
section IV to describe Hawking radiation. In section V we show how black hole entropy can 
be recovered from a microcanonical ensemble of states and discuss some issues raised by our 
approach in section VI. 



II. THE CLASSICAL LTB MODELS 

The LTB models describe self-gravitating dust. The energy- momentum tensor reads 
Tf_iu = e{T, p)Ufj,Uu, where = U^{t,p) is the four-velocity vector of a dust particle with 
proper time r and labeled by p (p thus labels the various shells that together form the dust 
cloud). The LTB line element is given by 

ds' = dr' - jf^M - R\p)dn\ (1) 
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Inserting this expression into the Einstein equations yields 
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where F{p) and E{p) are non- negative functions of the label coordinate and we have set 
SttG = 1. The case of collapse is described by drR{T,p) < 0. There still exists the freedom 
to rescale the shell index p, which we fix by demanding R{0,p) = p, so that for r = the 
label coordinate p is equal to the curvature radius R. Now we can express the functions 
F(p) and E{p) in terms of the energy density e{T,p) at r = 0. From (|), 



F{p) = / £(0,p)p2dp, 
Jo 

E{p) = hd,R{T = Q,p)f-- re{0,p)p'dp, 



(3) 



so that F{p)/2 may be interpreted as the active gravitating mass inside of the shell labeled 
by p and E{p) as the total energy inside the shell. An analysis of the classical solutions for 
these models can be found in 0. In the present work we discuss the canonical formalism 
for the so-called "marginally bound" models, defined by E{p) = 0. A generalization to 
non-marginally bound models is found in 0. 

A. Hamiltonian formalism 



Begin with the general ansatz for a spherically-symmetric line element, 



-N'dt' + {dr + N'dty + R'dn\ 



(4) 



where and iV^ are the lapse and shift function, respectively. The canonical momenta are 
given by 
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(5) 



where a dot denotes a derivative with respect to coordinate time t, while a prime denotes a 
derivative with respect to r. All variables are functions of t and r. A Legendre transformation 
from the Einstein-Hilbert action then leads to 



EH 



dt / dr ( PlL + PrR - NH3 - N'H^ ) + Sqi:, 



(6) 



in which the Hamiltonian and the diffeomorphism (momentum) constraint are given by 

32 ~ 



\ R 2Ry G 
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(7) 
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respectively, and the boundary action Sgs is discussed below. 

The total action is the sum of and an action S'^ describing the dust. The canonical 
formalism for the latter was developed in p], (see also ||^). The dust action reads 

/roo 
dt J dr {P,f - NH'^ - N'H^) , (8) 

where 



H^ = Pr\ll + ^ and H^ = r'Pr (9) 



are the dust Hamiltonian and momentum constraints respectively. 



B. Mass function in terms of canonical variables 



In the following we shall write the mass function F{p), which was introduced in (|^), in 
terms of the canonical data. This is essential for deriving consistent falloff conditions that 
are appropriate for a realistic collapse model. We begin by requiring the spacetime described 
by the metric (^) to be embedded in a LTB spacetime. Considering the LTB metric (0), a 
foliation described by functions r(r, t) and p{r,t) leads to 



N = 



R'^pp' - tt' 
^{fp'-pr'). 



(10) 



When these expressions for lapse function and shift vector into the expression for the canon- 
ical momentum Pi in (BI) and the equations (0) are used, we find 



R 



R'VT^ - Tt\ 



where T = \ — F j R. Solving for r' gives 



r 



and inserting this expression into ([T0|) yields 



T - 



LPr 



R 



which determines T according to 
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We can thus express F locally in terms of the canonical data as follows: 

" Pi R'^' 



(11) 



(12) 



(13) 



(14) 



(15) 
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This expression, although obtained here for marginal models, possesses a wider range of 
applicability, holding, in fact, for all cases 0. Further, it turns out that the functions, Pp, 
defined by 

and the mass function, F, form a conjugate pair of variables. Since R = F ai the horizon, 
J-" = there. We can check that though J-" appears in the denominator of (]T2|), r' is well 
behaved at the horizon: 



(R' + L) . 



(17) 



This is as it should be. We now make a canonical transformation in order to elevate the mass 
function F to a canonical coordinate. The canonical transformation, (r, R, L, P^, Pr, P^) — > 
(r, R, F, Pr, Pr, Pp), is generated by 



dr 



LPl - -RR' In 
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The action in the new canonical variables then reads 



EH 



dt / dr[ Prf + PrR + PpF -NH- N'Hr + S, 



(18) 

(19) 
(20) 

(21) 



where the new constraints are 



2L \ g:f 

H, = r'Pr + R'PR + F'Pp. 
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(22) 



We shall now discuss the boundary action Sqy. in niore detail. 



C. Boundary action 

Boundary terms are obtained from a careful discussion of the falloff conditions for the 
canonical variables, which were investigated in detail in [^. It turns out that the only 
boundary term is obtained from the variation of the hypersurface action with respect to L 
and reads 

y" dt A^+(t)5M+(t), (23) 

where iV+(t) = N{t,r — >■ oo) is the lapse function at infinity and M_|_(t) = F{r — )■ oo)/2 is 
the ADM mass. To avoid the conclusion that N^{t) is constrained to vanish, which would 
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freeze the evolution at infinity, the boundary term has to be canceled by an appropriate 
boundary action. This can be achieved by adding the surface action 

Sd^ = - j dtN+{t)M+{t). (24) 

Since varying A^_|_ would lead to a zero ADM mass, Kuchaf has argued in that A^_|_ has 
to be treated as a prescribed function. The lapse function gives the ratio of proper time to 
coordinate time in the direction normal to the foliation. Since N''"{r) vanishes for r — )■ oo, 
the time evolution at infinity is generated along the world lines of observers with r = const. 
If we introduce the proper time, f+, of these observers as a new variable, we can express the 
lapse function in the form N+(t) = r+(t). This leads to 

S9^ = - J dtM+f+. (25) 

and thus we have removed the necessity of fixing the lapse function at infinity. (In this 
is called 'parametrization at infinities'.) 

Out aim is to cast the homogeneous part of the action into Liouville form and to find 
a transformation to new canonical variables that absorb the boundary terms. This can be 
done by introducing the mass density F = F' as a new canonical variable and using the 
boundary condition F{0) = (which is appropriate for a collapse situation). Part of the 
Liouville form can then be rewritten as follows: 



9 = / dr PpSF - M+5f+ 
Jo 

= dr5T (y + ^ dr'P^(r')^ - 5(M+f+) (26) 

(see 0] for details). From (|2B|) we see that Pr = r+/2 + J^drPp. Thus Pr(oo) = f+/2. 
Thus the new action reads 

^EH = / ^ {^rf + PrR + Prt - NH^ - N'-Rs^ . (27) 

The constraints in the new variables are 



H, = r'Pr + R'PR-TP^. (28) 

For the Schwarzschild black hole, 2Pr is equal to the Killing time, T, in the exterior, as shown 
in [^. The Hamiltonian constraint can be greatly simplified if the momentum constraint 
is used to eliminate Pp = — Pf (see Appendix A of 0). The constraints (^) can then be 
replaced by the following equivalent set, 

H = (P,2 + j-i32) - ^ 0, (29) 
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Hr = r'Pr + R'Pr - TPr ^ 0. (30) 

These equations will be used as the starting point for the quantization in Sec. III. Another 
useful relation follows from ([12D , 



r = T ± y dR^^^y^ (31) 

and relates the dust proper time, r, to the Killing time, T, in the exterior of the collapsing 
dust ball, where the mass function is constant. The positive sign refers to contracting clouds 
and the negative sign to expanding clouds. 



D. Hamiltonian equations of motion 

Here we shall give the Hamilton equations of motion and derive Einstein's equation (^ 
from them. The Hamiltonian equations are generally given by 

X = {x,n[N] + nr[x']}, 

Px = {PxMN]+V-rW]], (32) 
where we have introduced the smeared constraints 

noo POO 

n[N] = / drN{r)H{r), H,[A^n = / drN''{r)Hrir). (33) 
Jo Jo 



Starting from the action (27), the Hamiltonian equations of motion are 



f = 2NPr + N''T', 
Pr = (NrPr)', 

R = 2NJ^Pr + NrR', 



r = (A^T)' 



= N— + iv^p' + r dm{f-) ( + — \ . (34) 

Consider now the momentum constraint in the form 

r.?^-'f^, (35) 
and use the Hamiltonian constraint 



— P / r2 



Note that 6F{r)/5T{f) = e{r - f). 
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so that 



Now using 

(0) can be compared to (O) and shows that 



r p 

= -Pp = (38) 



so that 



Pr = ^ (39) 



f 2NP^ r ' ^ ^ 



which we can solve for Pr. Inserting this expression into the Hamiltonian constraint 
gives 

4 4J^ 4J- ^ ' 

and solving for (drR)'^ leads to Einstein's equation (Q). Note that we did not have to specify 
the lapse function. 

The algebra of the constraints is not of the standard form (given, for example, in 0), be- 
cause we have used the momentum constraint to eliminate Pp in the Hamiltonian constraint. 
In fact, a short calculation gives 

{'H[N],'H[M]}=0, (42) 

{Ur [N^^] ,n[N]} = n [NrN'- - AT at;] , (43) 

{-HrWl -HrlM']} = Ur [[N\ M^] . (44) 

The Poisson bracket of the Hamiltonian with itself vanishes, in contrast with the general 
case in which it closes on the momentum constraint. The other brackets coincide with the 
general case. The transformations generated by the Hamiltonian constraint can thus no 
longer be interpreted as hypersurface deformations. They are in general not orthogonal to 
the hypersurfaces, but act along the flow lines of dust. 



III. QUANTIZATION 

We shall now apply the quantization procedure proposed by Dirac and turn the classical 
constraints into quantum operators, cf. Q. We begin with the expressions in (|29|) . Poisson 
brackets are translated into commutators in the Schrodinger representation by substituting 

Pr{r)-^-j{-y Pnir)^-j^y ^r(r)^--^ (45) 
I or(r) t oR[r) i ol (r) 



and having them act on wave functionals. The Hamiltonian constraint (^Oj) then leads to 
the WDW equation. 



(5r(r)2 5i?(r)2 ^ ' ^ ^ ' SR{r) 
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+ B{R,F)6{0r)- 



4J^ 



(46) 



where A and B are smooth functions of R and F that encapsulate the factor ordering am- 
biguities. We have introduced 5(0) in order to indicate that the factor ordering problem is 
unsolved and can be dealt with only after some suitable regularization has been performed, 
cf. |10]. That is, we choose the terms proportional to S{0) in such a way that the constraint 



algebra closes, which is usually called 'Dirac consistency'. Quantizing the momentum con- 
straint in (E^) by using (H3) gives 



+ R'- 



S 



r 
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^[t,R,T] = 0. 



(5r(r) ' " 6R{r) \6T{r) 
The next subsection is devoted to the application of a lattice regularization. 



(47) 



A. Lattice regularization 



For solutions of the constraints, we make the ansatz 



exp 



j dr r(r) >V(r,i?,F) 



(48) 



where VV(r, i?, F) is some function to be determined. It automatically satisfies the diffeo- 
morphism constraint. The Wheeler-DeWitt equation is second order in time derivatives so 
both positive and negative energy solutions exist, but we will confine our attention to the 
positive energy solutions above. It is worth noting that any functional 



vl/[r, R,T] = U (^-'- j dr T{r) W(r, R, F) 



(49) 



would satisfy the diffeomorphism constraint provided that W has no explicit dependence on 
the label coordinate r except through the mass function, F{r). We have chosen U = exp 
so that the wave-functional may also be factorizable on a spatial lattice, whose cell size we 
call a, taking cr — )■ in the continuum limit. Diffeomorphism invariance requires that the 
continuum wave-functional and all physical results be independent of the cell size. On the 
lattice, the argument of the exponential function becomes 0, 11 



j dr r(r) W{t,R,F) ^ a r,>V(r„ i?„ F,), 



(50) 



where Tj = T{rj), etc. This turns the wave-functional into a product state. 



nT,R,T] = l[i;,{T^,Rj,F^) = JJV^f exp 



(51) 



provided that U = exp. 
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Before proceeding further it is necessary to define what is meant by a functional derivative 
when functions are defined on a lattice The defining equations can be understood by 
analogy with the simplest properties of functional derivatives of the functions J{x) 



6j{x) 



5J{x) 



dyj{y) = 1 (52) 



and from these definitions follows 

6 



dyJ{y)m = 0(x). (53) 



5J{x 

On a lattice we define, for the lattice intervals Xi and Xj, 

|:^ = A(x.-x,) = lim^, (54) 

where Vi labels the i^^ lattice site and Sij is the Kronecker S, equal to zero when the lattice 
sites Xi and Xj are different and one when they are the same. Just as 6{y — x) is only defined 
as an integrand in an integral, so A{xi — Xj) should also be considered defined only as a 
summand in a sum over lattice sites. Hence 



lim -^a V J(r,) = lim a V = 1 (55) 

jjj^(^^Jiji)(i){ri) = lima^ A(ri -rj)0(ri) = 0(rj). (56) 

— - — --^-hm-— , (57) 
where Jj = J{xj). This is compatible with the formal (continuum) definition of the func- 



and 



It follows that 



tional derivative. 



B. Collapse Wave Functionals 



When (^J) and (|3) are applied to the Wheeler-DeWitt equation in ( ^61) and \E'[r, R, T] is 
taken to be a product state, one obtains an equation describing the wave functions at each 
lattice point 



drf 



+ A 



^dR, 



ipj ^ 0, 



(58) 



but there is a further restriction arising from the diffeomorphism constraint. Inserting the 
ansatz in (0) into ([79|) , we find 



OR, 



1 

J 



10 



err, 



92 w, 



5, 



0, 



(59) 



which must be satisfied independently of a. This is only possible if the following three 
equations are simultaneously satisfied at each lattice site M, 



dR 



drf 
B, = 0. 



dR] 



] _ 
'J J ^3 

3' 



OR, 



(60) 



Moreover, it is straightforward that the Hamiltonian constraint is Hermitean if and only if 



Aj = J^jdR^ In(mjlJ'jl), 



(61) 



where rrij is the Hilbert space measure. 

Unique solutions to the equations in (|60|) and having the form given in (|5TD have been 
obtained in all, even the non-marginally bound, cases 0. For the marginally bound models 
the solution for the phase VV, in the exterior, i.e., for shells that lie outside the apparent 
horizon {Rj > Fj), is 



(±) 



— tanh ^ — 



Zj > 1, 



(62) 



where Zj = ^jRj/ Fj. The positive sign refers to ingoing waves, traveling toward the horizon 
and the negative sign to outgoing waves, as can be seen from the signature of the phase 
velocity, 

% = (63) 

keeping in mind that Zj > 1. The phase velocity approaches zero at the horizon. In the 
interior, i.e., for shells that lie inside the apparent horizon [Rj < Fj), the solution is 



(±) 



Tj ± 2Fj [zj — tanh zj'j , zj < 1, 



(64) 



but here the the positive sign refers to outgoing waves and the negative sign to ingoing 
waves, traveling toward the central singularity, again as determined by the phase velocity. 
Furthermore, as shown in Appendix B of the system in (|60D determines not only Wj but 
the Hilbert space measure, xxij, as well. For the marginal models under consideration, trij is 
regular everywhere and given by 

m^- = Zj (65) 



up to a constant scaling. 
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IV. HAWKING RADIATION 



In this section we will argue that the states described above yield Hawking radiation. Our 
first approach will closely parallel Hawking's original work [|T2|. First we need to introduce 



the concept of a black hole into the formalism above, which, following ||13[, we do by taking 
the mass function to be of the form 



F(r) = 2Me(r) + /(r) 



(66) 



where B(r) is the Heaviside function and /(r) represents a dust perturbation (/(r)/2M ^ 
1). It is easy to see that with this choice of mass function, the black hole state factors out 
in ( |5lD and the remaining state then assumes the same form with F replaced by 2M and F 
replaced by /'(r). 



V^[r,i?,F] = e±^*^^« X exp 



±1 j drf'{r)Wf{T,R,M) 



def 



(67) 



where Wo = W(r(0), i?(0), F(0)) and the first exponent represents the black hole at the 
origin and the second, up to order /(r), represents the matter distribution that propagates 
in this background if we take F{r) ^ 2M in . 

Next, we must identify those quantum states that correspond to the ingoing and outgoing 
modes, respectively and evaluate an appropriate inner product. Since the description should 
refer to observers at infinity, the inner product will be evaluated on hypersurfaces of constant 
Killing (Schwarzschild) time T, and not on hypersurfaces of constant dust time (which 
corresponds to freely falling observers). It is not difficult to show that for contracting clouds 
the an infalling wave is given by 



/ = exp 



i j drf{r) (t + 8M (z - tanh"^ ^ 



being approximately 



^ exp 



i / drf{r) (T + %Mz) 



when T — )■ — oo and z oo and 

\E'^ ^ exp 



-i j drf'{r) ^T-SMtanh"^^ 



(68) 



(69) 



(70) 



when T — > +oo and z 1. Thus, the simple looking phase on Q"" scatters through the 
geometry to turn into the complicated looking phase on 53"'" near the horizon. This is similar 
to what happens in the geometric optics approximation. 

Equation (|69D represents infalling waves. We can think of it as a product over plane 
waves, one at each label r, as follows: 



n 



-iLjj[Tj+8Mzj 



(71) 
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This should represent a complete set of infalling modes at each label j, if we think of the Uj 
as the frequency of the modes. A complete set of outgoing modes on Q""^ would likewise be 
given by the functional 

j 



We now ask: what is the projection of our solution (^) on the negative frequency modes 
of the outgoing basis on Q""*". For this purpose we must consider the inner product of states 
on a hj^ersurface of constant Schwarzschild time T. It turns out that Hawking's thermal 
radiation is recovered if take the metric in the (r, R) plane is given by the quadratic term 
in the Hamiltonian constraint (^) instead of (|65D . We speculate that this has to do with 
the essentially classical role played by the black hole, described by ^'bh- Transforming to 
the metric in (i?, T) coordinates we find 



9rr 



R 



R-2M 



The required Bogoliubov coefficient is then given by the following inner product on a constant 
T hypersurface, 

Pif,u) = (v&Jv[;+) = f ^ dR, v[/+ V1/+, (73) 

which represents the negative frequency modes present in (^). A straightforward compu- 
tation then yields 



(in|iVo.t|m) = |/3(/,a;)p^n^ 



(74) 



which is interpreted as the eternal black hole being in equilibrium with a thermal bath at 
the Hawking temperature (SvrM)"^. Thus we have a functional Schrodinger picture for dust 
Hawking radiation (see fl^ for a generalization to the non- marginal models). 



An alternative approach, one that is better adapted to quantum collapse^ was considered 

in 



jT5]. By matching the shell wave functions describing gravitational collapse across the 
apparent horizon, it was shown that an ingoing wave on one side of the apparent horizon is 
necessarily accompanied by an outgoing wave on the other side. Furthermore, the relative 
amplitude of the outgoing wave is suppressed by the square root of the Boltzmann factor at 
the "Hawking" temperature of the shell. Strictly speaking the Hawking temperature, Th = 
(SttGM)"^, refers to an eternal black hole of mass M. The temperature appearing in the 
Boltzmann factor from matching shell wave functions across the horizon is Th = {4:7rGF)~^ , 
where F is the mass function and represents twice the mass contained within the shell. 
Diffeomorphism invariant wave functionals describing the collapse can also be matched and 
yields the same picture, but now the relative amplitude of the outgoing wave functional to 
the ingoing one is given by e~^^'^, where S is the entropy of the final state black hole. 
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V. BLACK HOLE ENTROPY 



As mentioned in the previous section, black holes with ADM mass parameter M are 
special cases of the solution in (|1]), obtained when the mass function is constant, F = 
2 GM, and the energy function is vanishing. This can be shown directly by a coordinate 
transformation of from the comoving system (r, p) to static coordinates (T, R), in which 
the metric has the standard from, 

ds^ = T{R)dT^ - J^-\R)dR^ - R^dQ\ (75) 

We imagine therefore that the eternal black hole is a single shell and represented by the 
mass function 

F{r) = 2Me(r) (76) 

(G = 1), where M is the mass at label r = and 0(r) is the Heaviside function. The mass 
density function is therefore 

r(r) = 2M5{r) (77) 

and, because of the 5— distributional mass density, the wave-functional in (^) turns into 
the waye-function 

^[t, R, r] = /o°° 'irr(r)W(r(r),i?(r),F(r)) ^ ^-iMWo{t,R,F) ^yg^ 

where r = t(0), R = R{0) and F = F{0). The Wheeler-DeWitt equation now becomes 
the Klein-Gordon equation describing the shell. Taking into account the factor ordering 
ambiguities and absorbing the M dependent term, which now renormalizes the potential, 
into the function B{R, F) we have 



+ J^TTT^ + A— + B 



dr^ dR^ OR 



e 



-iMWo{T,R,F) ^ (jQ^ 



In contrast with the case in which the mass density is a smooth function over some set of 
non-zero measure, no regularization is necessary here. This means that no further condi- 
tions must be met and therefore that the measure as well as the functions A{R, F) and 
B{R, F) will remain undetermined although the function A{R, F) will continue to be re- 
lated to the measure according to (^). Thus two conditions are required to proceed with 
the quantization of the black holes as described above. 

The first condition we impose is one on the measure appropriate to the Hilbert space of 
wave-functions. In the previous section we obtained the Hawking evaporation of a collapsing 
dust cloud surrounding a pre-existing black hole by taking the dust as a small perturbation to 



the black hole mass function in ([761). The calculation proceeded by evaluating the Bogoliubov 
coefficient in the near horizon limit outside the horizon and crucial to obtaining the correct 
Hawking temperature is the choice of measure appropriate for eternal black holes. The 
measure was obtained from the DeWitt supermetric, 'jab, on the configuration space (r, R) 
and can be read directly from the Hamiltonian constraint 
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It gives fi = 1/ v^l J^l, i.e., 

(vl/i,vl/2)= / -^^\^2 (81) 




as well as the function A{R, F) via the hermiticity condition (|65D. As long as J-" 7^ the 
Wheeler DeWitt equation can now be written as 



dr^ dm 



^ = 0, (82) 



where the positive sign in the above equation refers to the exterior, while the negative sign 
refers to the interior and i?^, is defined by 



R, = ± (83) 




The second condition arises because we are describing a single shell in this simple quantum 
mechanical model of an eternal black hole and because B{R, F) represents an interaction of 
the shell with itself. We simply demand there are no self interactions, i.e., that B{R, F) = 0. 
The quantum evolution is then described by the free wave equation in the interior, but by 
an elliptic equation in the exterior. This signature change has been noted in other models 
||16| , p!7| and occurs because of the behavior of J-", which passes from positive outside the 
horizon to negative inside. For the black hole, it means that its wave function is supported 
in its interior. The spectrum will be determined by the proper radius, Lh, of the horizon, 

dR 

Lh{M) = / (84) 



where Rh is its area radius. If we extend the coordinate R^ to range over (—00, 00), thereby 
avoiding any issues related to a boundary at the center, this simple model of a quantum black 
hole effectively describes a dust shell in a "box" of radius 2Lh{M), which itself depends on 
its total ADM mass. The stationary states describe a spectrum of the form (reintroducing 
h and G) 

AGM,Lh,j = Api (^j + , (85) 

where j is a whole number and Ap\ = hO is the Planck area. It is straightforward to show 
that Lh = TcGM, so that 

^ = ^Pi (j + I) , (86) 



where Aj is the horizon area |jT8| . 

As we show below, the equispaced area spectrum predicted by our simplified model of 
a quantum black hole implies that the entropy obeys the Bekenstein-Hawking area law 
provided that the area quanta are assumed distinguishable. The entropy therefore also 
admits a discrete and equispaced spectrum. What is the origin of the degeneracy that leads 
to the black hole entropy? Although the black hole is treated as a single shell with the 
spectrum in (^51), this single shell is in fact the end state of many shells that have collapsed 
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to form the black hole. Regardless of their history, we assume that each of the shells then 
occupies only the levels of (|85D , contributing some multiple of the Planck area to the total 
horizon area of the final state. A black hole microstate is thus a particular distribution of 
collapsed shells among the available levels. If the distribution of shells is such that Mj shells 
occupy level j, the black hole's total horizon area becomes 

J = ^piE(^ + ^)-^^- (87) 

and the (single shell) solution in (^) is to be interpreted as an excitation by A/" = Xlj -^j 
collapsed shells. 

The spectrum in ( ]87| ) represents the "area ensemble" and the number of black hole 
microstates giving the "area" A will depend on assumptions concerning the degeneracy 
of the microstates. Assuming the shells to be distinguishable, the number of states can be 
written in terms of the total number of area quanta, Q, and the total number of shells. A/", 

where 

Holding A fixed and extremizing the microcanonical entropy, Smicro = Ajslnfi, with respect 
to the number of shells gives 

Smicro = (2A;ij coth^i 75) , (90) 

which is in excellent (better than 96%) agreement with the Bekenstein-Hawking entropy. In 
addition to the exponential growth in the number of states, the area quantization in ( pT]) 
ensures that the entropy is effectively quantized in units of the Planck area, as originally 
proposed by Bekenstein . 

Note that it is quite a simple matter to show that had the shells been assumed indistin- 
guishable then the entropy would depend on the square- root of the area [Q. The fact that 
the area degrees of freedom must be treated as distinguishable runs contrary to our intuition 
for elementary degrees of freedom in quantum field theory and calls into question whether 
"area" is a fundamental quantity in quantum gravity. 



VI. DISCUSSION 

In our contribution we have described a canonical quantization of collapsing, inhomo- 
geneous dust and some of its implications. While here we have confined ourselves to four 
dimensions and a vanishing cosmological constant, our set up can be extended to describe 
dust collapse in any dimension both with and without a cosmological constant and in all 
cases it is possible to show that static black holes will radiate at the appropriate Hawking 



temperature of the black hole [21, 021, at least in the semi-classical approximation in which 
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the mass function is taken to be of the form given in (^). Two, hkely related, features of 
this description remain to be understood. Firstly, why does the Hawking picture rely on the 
measure derived from the classical Hamiltonian constraint to describe eternal black holes? 
Secondly how does the semiclassical picture hold up during collapse, i.e., when the dust is 
not taken to be a perturbation on the background of a massive black hole? In principle, the 
second question could be addressed by pursuing the approach of [|T3|, but it still remains to 



be done. Another related problem concerns the regularization used to define the functional 
Wheeler-DeWitt equation (|46|) . The lattice regularization we have used is an ad hoc regu- 
larization in which the divergent terms have to cancel each other. It leads directly to (|60D, 
implying that the Hamilton Jacobi solution is exact and, if factorizability on the lattice is 
required, making (|62|) and (|6^ unique as well. However, regularization independence of the 
results described here remains an open question and factorizability on a spatial lattice may 
be too strong a condition to impose. 

Fundamental questions concerning the nature of the quantum gravitational degrees free- 
dom remain, and these are most manifest in the assumptions that must be made concerning 
the statistics obeyed by them. For the BTZ black hole in 2+1 dimensions it is necessary 
to employ Bose statistics instead of Boltzman statistics in order to recover the entropy p3 |. 



This corresponds to the fact that the BTZ black hole admits an equispaced mass spectrum 
(as opposed to an area spectrum) and must be treated in the energy ensemble. Its heat 
capacity is positive, whereas it is negative for the Schwarzschild black hole. The actual 



computation of the entropy is similar to that employing the AdS / CFT correspondence |24 
and the entropy is found to depend on two quantities, viz., the energy of what is taken to be 
the vacuum solution, Aq, and a constant, Mq, arising out of a boundary contribution at the 
origin.^ A comparison between the result from the canonical theory and the result obtained 
via the AdS / CFT correspondence yields an effective central charge 



1 



(91) 



which must be set to 3l/2Gh to achieve agreement with the Bekenstein-Hawking entropy. 
This is the central charge of the Liouville theory induced at spatial infinity by 2+1 dimen- 



sional gravity 



In three or more spatial dimensions, with a negative cosmo logical constant, the spectrum 



of states describing eternal black holes has a more complicated description |2^. For small 
black holes, by which we mean black holes whose horizon radii are much smaller than the 
AdS length, an approximate area spectrum, similar to (p6|), is obtained and the entropy 
is recovered by working in an area ensemble in which the degrees of freedom are assumed 
distinguishable. In this case, the black hole heat capacity is negative. On the other hand, 
for large AdS black holes, i.e., black holes whose horizon radii are much larger than the 



Whereas a boundary contribution from the origin in 3+1 dimensional collapse would represent a singular 
initial configuration and is therefore set to zero, in 2+1 dimensions a non-vanishing contribution is essential 
to allow for an initial velocity profile that vanishes there. This does not lead to singular initial data and 
the boundary contribution does not have the interpretation of a point mass situated at the center . 
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AdS length, it is the mass that is quantized in integer units and the (mass) spectrum turns 
out to be independent of the gravitational constant, G. This is similar to what happens 
with the BTZ black hole (for which, however, there is no "small black hole" limit). The 
entropy is now obtained in an energy ensemble and the degrees of freedom must be assumed 
indistinguishable. It can then be shown that the thermodynamics in the large black hole 
limit is inextricably connected with the thermodynamics in the opposite limit by a duality 
of the Bose partition function. The gravitational constant, G, absent in the mass spectrum, 
reemerges in the thermodynamic description via this duality and the black hole heat capacity 
is positive. It appears that the Hawking-Page transition separates the particle-like 
degrees of freedom of large black holes, which must be counted as indistinguishable, from 
the geometric degrees of freedom of small black holes, which are counted as distinguishable. 
This remains rather mysterious and seems worth pursuing further, but that is complicated 
by the fact that it is difficult to obtain closed form solutions for the spectrum of AdS black 
holes between the two limits described. 

In conclusion, we now seem to be in a position to address some of the issues that have 
surrounded Black Hole thermodynamics for decades, such as the information loss problem, 
questions about singularity avoidance and the radiation from a naked singularity, at least 
within the context of a special class of models. We expect progress in addressing these and 
the new issues that have arisen since developing the quantization described here to be made 
in the near future. 
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